Role of quantum coherence in chromophoric energy transport 
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The role of quantum coherence and the environment in the dynamics of excitation energy transfer is not fully 
understood. In this work, we introduce the concept of dynamical contributions of various physical processes to 
the energy transfer efficiency. We develop two complementary approaches, based on a Green's function method 
and energy transfer susceptibilities, and quantify the importance of the Hamiltonian evolution, phonon-induced 
decoherence, and spatial relaxation pathways. We investigate the Fenna-Matthews-Olson protein complex, 
where we find a contribution of coherent dynamics of about 10% and of relaxation of 80%. 
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I. INTRODUCTION 

Exciton transfer among chlorophyll molecules is the energy 
transport mechanism of the initial step of the photosynthetic 
process. Light is captured by an antenna complex and the ex- 
citon is subsequently transferred to a reaction center where 
bio-chemical energy storage is initiated by a charge separa- 
tion event JD. This transfer process has been studied us- 
ing classical Forster theory or a (modified) Redfield/Lindblad 
description El[2ll|5jSll3[8l|E[l0]- Several measures 
such as the energy transfer efficiency/quantum yield, trans- 
fer time, and exciton lifetime have been employed to eluci- 
date the performance of exciton transfer ifTTl [T2l [T3ll . Re- 
cent experiments suggest evidence of long-lived quantum co- 
herence in the Fenna-Matthews-Olson (FMO) protein com- 
plex of the Green-Sulphur bacterium Chlorobium tepidum and 
in the reaction center of the purple bacterium Rhodobacter 
sphaeroides lH4l[T5l . A naturally arising question is the role 
of coherence in the biological function of the aforementioned 
chromophoric complexes. 

In this work, we investigate relevant quantum coherence ef- 
fects by an in-situ analysis of a success criterion for the initial 
step in photosynthesis, the energy transfer efficiency (ETE). 
The dynamics of an excitation in multi-chromophoric com- 
plexes can be described in terms of an environment-assisted 
quantum walk [ 10 1 : 

±p(t)=M P (t). (1) 

The evolution generated by the superoperator M. connects the 
population and coherence elements of the density matrix p. 
We assume in this paper that M. is time-independent. We 
would like to explore and characterize the dynamics of Eq. ([IJ, 
especially the role of quantum coherence, the environment, 
and spatial energy transfer pathways. For closely packed mul- 
tichromophoric arrays, such as the FMO complex, one has to 
account for strong inter-molecular coupling and quantum co- 
herence effects. A classical approximation of the master equa- 
tion would be an insufficient description of the open quantum 
dynamics. Specifically, we want to avoid any comparison of 
the actual open quantum system under study with a fictitious 
or abstract model system, such as a high temperature limit, a 
strong decoherence model, or a semi-classical Forster method. 



This is in contrast to studies e.g. in the area of quantum infor- 
mation that compare the quantum dynamics to classical dy- 
namics, for example in the case of the comparison of a classi- 
cal random walk to a quantum walk fl6l . 

We quantify the role of the various physical processes in- 
volved in the energy transfer process in terms of their con- 
tribution to the ETE. Formally, we will partition the overall 
ETE, 77, into a sum of terms, 

V = ^2vk, (2) 

k 

corresponding to a physical decomposition A4 = ^2k-M.k- 
Each term rjk can be interpreted as a contribution to the over- 
all efficiency originating from a particular process Mk- For 
example, we will split the superoperator into the major com- 
ponents describing the exciton dynamics: coherent evolution 
with the excitonic Hamiltonian, relaxation within the single- 
exciton manifold, and dephasing. The r\% associated with the 
coherent part will then give an indication of the role of quan- 
tum evolution to the energy transfer efficiency and hence to 
the biological function within a particular chromophoric com- 
plex. We note that the exact partitioning of the ETE into a sum 
of terms like Eq. Q is a non-trivial task: as will be described 
below, the ETE essentially involves an exponential mapping 
of the complete superoperator M.. A separation of the ETE 
into a product of terms would seem more natural but would 
not allow the interpretation of t]k as contributions. 

In the following sections, we briefly discuss the structure of 
the superoperator M. and introduce two complementary mea- 
sures of efficiency contributions: one is based on a Green's 
function method and the other is derived from energy transfer 
susceptibilities. We apply these two approaches to the study of 
the ETE in the FMO complex. We employ a standard Redfield 
model with the secular approximation which leads to a mas- 
ter equation in Lindblad form |[T9l . This model captures major 
decoherence effects such as relaxation and dephasing. We also 
include spatial correlation of the fluctuations. The Markovian 
approximation neglects temporal correlations in the phonon 
bath which can be relevant in photosynthetic systems and will 
be treated in subsequent work. We believe that this model 
and our methods can provide insight into the role of quantum 
coherence in energy transfer, a process that occurs in noisy, 
ambient temperature environments. 
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II. MASTER EQUATION FOR MULTICHROMOPHORIC 

SYSTEMS 

The transport dynamics of a single excitation is described 
by a master equation for the density matrix that includes co- 
herent evolution, relaxation, and dephasing. Moreover, the 
exciton can recombine or be trapped. The Hamiltonian for an 
interacting TV-chromophoric system in the presence of a sin- 
gle excitation can be written as |8 1: 



iV 



iV 



Hs=J2 € ™\™>){m\ + V mn(\m)(n\ + \n)(m\), (3) 

m— 1 n<m 

where the Hilbert space basis states \m) denote the presence 
of an excitation at the mth chromophore and e m are relative 
site energies with respect to the chromophore with the low- 
est absorption energy. The V mn can be due to Coulomb cou- 
pling of the transition densities (Forster) or due to overlap 
of electronic wave functions (Dexter). We denote the eigen- 
basis of the Hamiltonian ([3} as the exciton basis \M) — 
J2 m c m(M)\m), where HgjM) — £m\M). The multichro- 
mophoric system interacts with a thermal phonon bath. The 
dominant component of the system-bath Hamiltonian is as- 
sociated with site-energy fluctuations ifTTl Il8l . i.e. Hsb = 
^ m g m |m)(m|, where q m are operators describing the cou- 
pling to the coordinates of the harmonic-oscillator bath. The 
phonon terms q m \m)(m\ induce relaxation and dephasing 
without changing the number of excitations. We assume 
that the bath correlator can be simplified as (q m (t)q n (0)} = 
C m n(q(t)q(0)}. C mn is a dimensionless time-independent 
factor that takes into account the spatial correlations in the 
phonon bath. For spatially uncorrelated environments it will 
simply be given by C mn = 8 mn . In this work, we will also 
take into account a phenomenological model for these corre- 
lations as will be explained later. The time-dependent part of 
the correlator is the same for all sites IfTTl . Additionally, there 
are two processes that lead to irreversible loss of the exciton 
iflOl ITTl [T2l [TJft . One is the excitation loss due to recom- 
bination of the electron-hole pair. The other mechanism de- 
scribes the excitation transfer to the reaction center (acceptor) 
and subsequent trapping associated with the charge separa- 
tion event. These effects are taken into account by the anti- 
Hermitian Hamiltonians, — ii? rcC omb = — i^I 1 Ylm \m)(m\, 
with r the inverse lifetime of the exciton and — \H\ 



-i/i^ K m |m) (m|, with K m the trapping rates at site to. 

In summary, the dynamics of the reduced density matrix of 
the system can be described by the Lindblad master equation 
in the Born-Markov and secular approximations as 1 19 1: 



trap 



dp® 

dt 



--[H s + H LS ,p(t)] + Cp(t) (4) 

~^{H Iccomh ,p(t)} - ^{H trap ,p(t)}, 

where { , } denotes the anti-commutator. The right-hand side 
of Eq. (|4} defines the superoperator M.. £ is the Lindblad 
superoperator derived from the phonon bath coupling, 



£p(t) = ]T lmn (cu)[A m (cu)p(t)Ai(cu) 



(5) 



-~A m (cj)Al(u)p(t) - l -p{t)A m {uS)A\ (w)]. 

The sum runs over all possible transitions in the single 
exciton manifold and all the sites. The Lindblad gener- 
ators are A m (cu) = E eM - £N= ^ C Wc m (iV)|M)(iV|, 
where the summation runs over all transitions with fre- 
quency uj in the single-excitation manifold. The Fourier 
transform of the bath correlation function leads to the 
rates 7„ m (cj) = 2nC mn [J(uj)(l + n(uj))+ J(-uj)n(-uj)], 
where n(ui) is the bosonic distribution function at temper- 
ature T. Here, we assume an Ohmic spectral density with 
J(oj) = for u < and J(uj) = f^expf-^) else- 
where, with cutoff oj c , and reorganization energy i?R = 
H J °° duj^^p- 1171 . The pure dephasing part of Eq. (5b is 
obtained in the limit uj — > 0. The Lindblad generators are 
An(0) = Em \c m (M)\ 2 \M)(M\ and the rate 
2nC mn j^pkT. The Lamb shift Hamiltonian [19] is given by 

Hls = E R J2M, m \c m (M)\ 4 \M)(M\, where we only take 
into account the most significant, diagonal part, see Ref. lfT8l 
and supplementary material thereof. 

The competition of trapping and recombination processes 
leads to the concept of the energy transfer efficiency (ETE) as 
the integrated probability of successful exciton entrapment by 
the reaction center lfT0lflTl[T3l . 



77 = 



y o Ti-{H tiapP (t)}dt. 



(6) 



The efficiency is suppressed by the lifetime of the excitation 
due to exciton recombination. Quantitatively, the suppression 
is given by the "deficiency", fj = |/ °° Ti{H Iccomh p(t)}dt. 
One finds the relation r\ + fj = 1 which implies that the exci- 
tation ultimately either is trapped or recombines. 



III. CONTRIBUTIONS TO THE ENERGY TRANSFER 
EFFICIENCY 

We partition the efficiency based on a Green's function 
method. First, note that with Eq. <[TJ the efficiency can be 
simplified as 



V 



-Tr{ff trap .M-V(0)}, 



(7) 



The trapping and recombination Hamiltonians ensure that M. 
is invertible since non-trivial eigenstates with zero eigen- 
value, such as the usual thermal equilibrium, do not ex- 
ist. The Green's function interpretation of —Mr 1 can 
be readily confirmed by examining the Laplace-transformed 
master equation. Next, the superoperator is decomposed 
into M. = Href + 1Z, where the recombination-trapping part 
= -^{#rccomb,/o(£)} - H^ap , p(t) }, is taken as a 
reference point and the remainder is given by IZp = — i [Hs + 
HhSi p] + £p- This allows us to express the Green's function 
exactly as 



M- 



n 



ref 



(8) 
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The first term in ([8]) describes the irreversible dynamics re- 
sulting from excitonic transfer to the acceptor and recombina- 
tion of the excitation. The second term accounts for the effect 
of the Hamiltonian and the Lindblad operator. This term can 
be decomposed into various physical processes such as co- 
herent hopping or phonon-mediated jumps from site-to-site, 
7Z = ^2 k TZk- We use Eq. and the relation T] = 1 — f\ to 
arrive at the partitioning of the efficiency, rj — X)fc with 

Vk = -^Tr{ J ff rccomb H re ^ fc M-V(0)}. (9) 

The part related to the first term in Eq. <|8j gives just one 
and thus cancels with the one from the conservation property. 
Note that each r)k is dependent on the initial state, the over- 
all dynamics Mr 1 , and the particular process IZ^ weighted 
by a diagonal trapping-recombination operator H IecoIa bH~^. 
The trace then takes into account the local effect of IZk on the 
population of all the site basis states \m). 

IV. APPLICATION TO THE FMO COMPLEX 

The Fenna-Matthews-Olson protein complex is a chro- 
mophoric trimer where each subunit consists of seven chloro- 
phylls embedded in a protein environment ifPfl |2TI l22l . It 
connects the photosynthetic antenna with the reaction cen- 
ter in Chlorobium tepidum. We use the seven-level Hamil- 
tonian of the FMO complex as given in Ref. 1T71 and the 
bath spectrum described earlier with the reorganization en- 
ergy En = 35cm _1 and cutoff uj c — 150cm" 1 , inferred from 
Fig. 2 of Ref. |[T8l . The transfer of the excitation from the 
FMO to the reaction center occurs via site 3 with the rate k%, 
which is a free parameter in our simulations and, if not oth- 
erwise stated, is taken to be k 3 = lps -1 iflOl . The exciton 
lifetime is assumed to be l/T = Ins ||231 . The initial state 
for the simulation is taken to be an unbiased classical mixture 
of all sites except the trapping site 3. At first, we will assume 
that the phonon modes at different chlorophylls are not spa- 
tially correlated with each other, which is equivalent to setting 
Cmn — 5 mn . Later, we will relax this assumption and inves- 
tigate the contributions as a function of a phenomenological 
parameter, the correlation length R c . 

In Fig. [T] (left panels) we present the contributions to the 
ETE as a function of three relevant system parameters for the 
three fundamental physical processes: Hamiltonian evolution 
with Hs and exciton relaxation and dephasing, both contained 
in L. The parameters explored are the reorganization energy 
En, which is a linear prefactor of all the decoherence rates, 
the transfer rate to the reaction center k 3 , and temperature 
T. In Ref. Q~0l|, we reported the overall behavior of the ef- 
ficiency with respect to these parameters, especially the en- 
hancement of the ETE by about 20% compared to no phonon- 
bath, E-r = 0. Even more significantly, the transport time of 
the excitation to the reaction center reduces from over 50ps 
to under 5ps. Here, by providing a more detailed analysis 
of the underlying processes, we observe a crossover from a 
purely quantum regime to a relaxation-dominated regime with 
increasing phonon bath coupling strength En, see Fig.[T](a). 



At room temperature and En — 35cm _1 the purely quan- 
tum mechanical contribution of coherent hopping due to the 
Hamiltonian is around 10%. This is the main result of this 
work and is to be understood within the context of our model 
and its approximations. The major contribution is due to ex- 
citon relaxation induced by the phonon-bath coupling. Site 3 
has a large participation in the lowest-energy exciton: relax- 
ation helps the transport from the energetically higher initial 
state towards that site. Dephasing processes have a contribu- 
tion of around 8% (see explanation below). To summarize, 
the results so far indicate that the process of energy relaxation 
determines the high efficiency and the fast transfer times f the 
Fenna-Matthews-Olson protein complex and thus is essential 
to the biological function of this system. Supporting this state- 
ment is a recent study that finds that the protein environment 
induces a static red-shift of the energy of site 3 11221 . 

Slow movement of the protein scaffold will lead to static 
disorder in the Hamiltonian and thus to inhomogeneous 
broadening of spectral lines. We obtain the effect of static 
disorder on the efficiency and the contributions by averaging 
over diagonal disorder of the site energies e m . We assume 
that the site energies are normal distributed around the Hamil- 
tonian given in ifTTl and with a site-dependent FWHM a m . 
Based on the supplementary material of Ref. [22 1 we choose 
°"i,3,4 = 60cm~ 1 , (72 = 100cm -1 , and (75,6,7 = 120cm _1 . 
These values lead to good agreement of measured and simu- 
lated spectra and can be explained by the different amounts of 
mobile water molecules in the vicinity of the different chloro- 
phylls of the FMO protein complex. In Fig. [T]the resulting 
distributions of the efficiency and the contributions are de- 
picted as error bars. At small En the efficiency and the preva- 
lent quantum mechanical contribution have a broadening of 
around 10%. This is due to the fact that the static disorder 
changes the energy difference of interacting sites which can 
facilitate or inhibit transport. At larger En the efficiency has 
a small broadening of less than 1%. In the presence of an en- 
vironment the system is more robust against static disorder. A 
smaller contribution of the quantum mechanical part is com- 
pensated by a larger contribution of relaxation and vice versa. 

Fig.[T](b) shows the dependence of the contributions on the 
transfer rate to the acceptor. The behavior can be explained 
by the characteristic time scales of the various processes. At 
room temperature and En = 35cm _1 , the relaxation pro- 
cesses occur on a time scale of lps and are most important 
when the acceptor rate is of the same order of magnitude. The 
Hamiltonian and dephasing act on faster time scales such that 
their contributions are stronger at higher acceptor rate values. 
The contribution of dephasing can be understood from the fact 
that, within the Redfield/Lindblad treatment, dephasing acts 
in the energy basis: an initial state localized at a site dephases 
to a state that has overlap with the target site. For example, 
the initial state p(0) = |1)(1| = £m,n Cl (M)cl(N) \M) (N\ 
evolves in the presence of only dephasing to p(t — > 00) = 
J2 M \ci(M)\ 2 \M)(M\, which has (3|p(f -> 00) |3) ^ 0. 

Finally, Fig.[T](c) shows the temperature dependence of the 
ETE contributions at En = 35cm _1 . The role of the Hamil- 
tonian is about 20% at zero temperature. The relaxation pro- 
cesses are only weakly dependent on T: the main contribu- 
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FIG. 1: Contributions of the basic physical processes to the overall energy transfer efficiency of the Fenna-Mathews-Olson protein complex 
according to Eq. {5} (left panels) and Eq. l[TTJ with normalization \12\ (right panels). The physical processes are the Hamiltonian evolution 

(red, ), exciton relaxation (green, — ), and exciton dephasing (blue, ■ ■ ■)■ They add up to the efficiency (black, — ) (left) or 1 (right). 

Depicted are the contributions as a function of a) reorganization energy, essentially the overall strength of the phonon-bath coupling b) the 
transfer rate to the reaction center, and c) the temperature. We assume that the initial state is an unbiased classical mixture of all the sites 
of the FMO complex except target site 3. Typical parameters (shown as vertical lines) are En = 35cm - , T = 295K, = lps -1 , and 
r = Ins -1 [ 1 8 1 . Using these parameters, the first (second) measure on the left (right) gives an estimated contribution of 80% (87.5%) for 
exciton relaxation and 10% (7.5%) for the Hamiltonian evolution at room temperature. 



tion is the temperature-independent spontaneous emission of 
energy into the phonon bath, leading to energy funneling to- 
wards the site with the lowest energy, site 3. The behavior 
of the dephasing processes is explained by 7^ ~ kT. In all 
figures, the contribution of the Lamb shift Hamiltonian is less 
than 1.5% and thus not depicted. 

Up to this point we assumed that the site energy fluctuations 



of a particular site caused by the coupling to a vibrational bath 
are uncorrelated from the fluctuations at another site. We relax 
this assumption to take into account spatially correlated fluc- 
tuations present in realistic chromophoric systems embedded 
in a protein environment. Specifically, instead of an uncorre- 
lated environment with C mn — S mn we include correlations 
with C mn = e~ Rmn / Ro |[T8l . R mn is the intermolecular dis- 
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FIG. 2: Efficiency and contributions of the basis physical pro- 
cesses as a function of the correlation radius in the environment 
in the Fenna-Matthews-Olson complex. The parameters used are 
E R = 35cm" 1 , T = 295K, k 3 = lps~\ and T = Ins" 1 . The 
lower panel is a magnification of the physically relevant region of 
the upper panel. 



tance between chlorophyll m and n and R c is the correlation 
radius. In the limit of R c — > one has the uncorrelated case. 
In the limit of R c — > oo all fluctuations are perfectly corre- 
lated and thus lead only to an irrelevant global phase in the 
exciton dynamics: this limit is equivalent to no phonon bath 
at all. In Fig. [2] we plot the efficiency and the contributions as 
a function of the correlation radius for the Fenna-Matthews- 
Olson protein complex. We recover the two limits: at small 
correlation radius the efficiency and the contributions are the 
same as in Fig. [T] For large correlation radius the quantum 
coherence contribution is dominant and there is no contribu- 
tion of relaxation and dephasing. In between, we observe a 
crossover which happens when the correlation radius is of the 
order of the intermolecular distance. 



V. BEYOND LOCAL CONTRIBUTIONS 

The method for partitioning the efficiency into contribu- 
tions presented in the last section sheds light on the role of 
the basic physical processes in the overall open systems dy- 
namics. Our approach resembles the sojourn expansion used 
in Ref. [ 1 1 1 for a classical random walk and the method used 



in Ref. [12| for studying incoherent versus coherent energy 
transfer in the Haken-Strobl model. However, it shows a cer- 
tain local feature: only processes that are directly connected 
to site population elements play a role. This can be seen in 
Eq. ^ from the fact that Ti~ e ] and .ffrecomb are diagonal op- 
erators. Thus the contributions of all IZk which do not lead 
to transfer to a site population vanish in the trace, e.g., in the 
case of coherence to coherence processes. Note that this ef- 
fect would have been pronounced had we not used the identity 
rj = 1 — fj before Eq. |9]). In this case only those TZj. that cre- 
ate overlap with the target site 3 would contribute. To over- 
come these issues, we develop our complementary measure 
motivated by the concept of energy transfer susceptibilities, 
discussed in 1101 . 

In this method, we partition the efficiency into terms involv- 
ing the variational change of the excitation trapping probabil- 
ity density j^Tr {H trap p(t)} with respect to the different phys- 
ical processes. We first provide a representation for a quantum 
master equation in the Markov approximation that is tailored 
to the problem at hand: 



dt 



P(t) 



i r t 

V- / F(t,t')M k F{t',0)p(0)dt'. (10) 
i. 1 Jo 



We used the decomposition M. = Y^k ^ k an< ^ tne 1 uan ~ 
turn map p(t) = J-(t, 0)p(0). This equation is equivalent to 
Eq. ([lj, but the individual terms in the sum are different from 
just A4k- the integrand can be understood as the realization of 
a non-unitary quantum walk in Liouville space. The quantum 
walk starts at the initial state p(0), evolves with T for a time 
interval of [0, t'], then a perturbation M. k is applied at an ar- 
bitrary time t', after which it undergoes a post-evolution until 
time t. The second evolution !F(t, t') leads to quantum inter- 
ference between Ai k and other generators. Integrating over 
all possible paths gives the average effect of the kth generator 
to the time variation of the density matrix. 

Eq. ( 10 1 can be reexpressed by introducing scalar dimen- 



sionless quantities associated with each term M. k , such 
that J\4k — » \k-M-k m the neighborhood A^ — ► 1. This leads to 
the master equation Jj/?(t) — \ ^2 k gf^p(t), which is equiv- 



alent to Eq. ( 10 1. We employ it to partition the ETE in Eq. d6j) 



into contributions given by 



Vk 



dt 



dt '^d^ TT{Ht ^ 



Pit')}- 



(11) 



We assume that the initial state p(0) has no overlap with the 
trapping sites. The terms | -r^Tr {H trap p(t)} can be inter- 
preted as the susceptibility with respect to the process M. k 
of the exciton trapping probability density |Tr{_ff t rapP(^)} at 
time t. The double integration is then to be considered as a 
time averaging of this quantity. The factor j in the integrand 

I by numerical inte- 
are not normalized, 



arises from Eq. (IOi. We solve Eq. (11 
gration. The resulting contributions % 
meaning that < r\ k < 1 does not necessarily hold. To define 
a proper normalization we separate the r\ k into positive and 
negative terms, i.e. 77 = 77+ + 77" = J^k + Vk + + J2k- Vk- 

< and the overall posi- 
Now one can introduce a 



where one has 77^ > and 77^, 
tive or negative contributions 77 
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FIG. 3: The Fenna-Matthews-Olson protein complex and the contri- 
butions in percent of the relaxation pathways to the overall relaxation 
contribution of the energy transfer efficiency, using Eq.[TT]and the 
normalization Eq. [T2] For clarity some contributions below 2% are 
not depicted. The initial state is a classical mixture of all sites except 
target site 3. Standard parameters are Sr. = 35cm -1 , T = 295K, 
— lps -1 , and F = Ins -1 . The contributions show directionality 
in their sign and reveal the important pathways. 

normalization given by: 



This allows us to interpret rjk ± as positive or negative percent- 
age contributions to the overall positive or negative contribu- 
tion. The normalized contributions rjf. ± show similar behav- 
ior for the basic processes as the Green's function measure 
Eq. d9), see Fig. [T] (right panels). This can be seen as ev- 
idence that, one the one hand, the normalization procedure 
was appropriate and, on the other hand, that despite its lo- 
cal feature the Green's function measure obtains consistent 
results. The second measure generally assigns less contribu- 
tion to the free Hamiltonian, about 7.5% at En = 35cm _1 



and T — 295K. For these parameters, the contribution of 
relaxation is about 87.5% and dephasing about 5%. At zero 
temperature the Hamiltonian contribution is 12.5%. 

The measure based on susceptibilities can be used to quan- 
tify the contributions of exciton relaxation pathways in the site 
basis. Formally, we look at the contribution of all site-to-site 
jumps and the corresponding damping of diagonal popula- 
tions, see Ref. ifTOl . In Fig. [3] we show the contributions of 
the various pathways for the Fenna-Matthews-Olson protein 
complex when the system is initially in a classical mixture of 
all sites except the target site 3. Jumps toward the target site 
3 contribute positively while jumps away from the target site 
contribute negatively. Large contributions come from nearest 
neighbor jumps, and site 4 and 7 are revealed as hubs in the 
transfer toward site 3. 



VI. CONCLUSION 

In this work, we have addressed the role of quantum co- 
herence and the environment in excitonic energy transfer. To 
this end, we have characterized the underlying processes con- 
stituting the open quantum walk of the excitation in terms of 
their contribution to the transfer efficiency. The methods pre- 
sented here are general and can be applied to a large class of 
transport systems in the presence of Markovian environments. 
Within both the Green's function and the energy transfer sus- 
ceptibility formalisms we conclude that the major part of the 
high efficiency of the Fenna-Matthews-Olson protein com- 
plex of about 80% or 87.5% is due to environment-induced 
relaxation down to the lowest energy site. The role of quan- 
tum coherence induced by the Hamiltonian dynamics can be 
quantified at around 10% or 7.5% respectively. Furthermore, 
we used the susceptibility measure to assign percentage-wise 
contributions to exciton relaxation pathways in the molecular 
basis. The detailed analysis of the open quantum dynamics 
presented in this work could be harnessed for engineering ar- 
tificial materials such as quantum dots 11251 to achieve optimal 
energy transport in realistic environments. 
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